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Abstract 

We point out that both superexchange between Cu02 layers, and interlayer 
tunneling, derive from frustrated one-particle hopping between layers, and 
that they should be treated on an equal footing. Doing so, we arrive at 
a new view of the nature of pairing in the cuprate supercondcutors, which 
explains the striking even <-> odd selection rule observed by Keimer in neutron 
scattering by YBCO. 
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Recently, two different hypotheses have been put forward for the "spin gap" phenomenon 
exhibited by the bilayer structures of cuprate superconductors. (As in YBCO and BISCO 2- 
layer materials). Millis and Monien and co-workers |^ have proposed that what is occurring 
is pair formation between the two layers with one of the pair on each layer, motivated by 
the interlayer antiferromagnetic superexchange. Strong and myself have proposed that 
the interlayer pair tunneling Hamiltonian produces a correlated state with pairing on each 
layer for each momentum k, but the different fc's are independently phased. Both postulate 
"preformed BCS pair" states as the essential nature of the spin gap phenomenon. 

Although it is not explicitly mentioned in the former paper, both theories have one vital 
element in common: the coherent motion of single electrons between the two layers must be 
blocked. As I have emphasized elsewhere, antiferromagnetic "superexchange" interaction 
exists only if the single-particle kinetic energy is frustrated, and is a consequence of the 
second-order, virtual action of this kinetic energy, causing hopping of electrons between the 
relevant sites. If the kinetic energy were diagonalized, the Fermi surface would be split 
into separate surfaces for even and odd linear combinations of the two layers, which, being 
orthogonal, would exhibit only ferromagnetic exchange interactions. (The reader should note 
that although the rest of this article will be couched in language referring only to bilayer 
materials, only a slight modification allows one to present similar arguments for coupled 
layers either single or multiple layers.) 

In fact, the Millis-Monien et al papers miss another feature of interlayer superexchange. 
The hopping matrix element t± between layers is of course diagonal in k\\, the interlayer 
momentum, so that superexchange, ~ t^/U, must also contain an extra momentum S- 
function. That is, 

-ly _ \^ +(1) J2) +(2) (1) /..v 
k k' 

xtk ■ t^/ /energy denom 

and the only way to satisfy the pairing condition is for k' = —k. Thus these two mechanisms 
share the characteristic fc-diagonality, which leads to approximate independence of different 
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parts of the Fermi surface; and the Strong-Anderson mechanism for producing a spin gap is 
common to both. 

A third relevant theory paper is the weak-couphng solution of the two-chain Hubbard 
model by Balents and Fisher Q], in which much of the relevant region of parameter space 
seems to exhibit the same fixed point with a spin gap and one gapless charged excitation, 
which seems to interpolate between the pairing schemes of Ref. (1) and (2). 

We argue here that the two mechanisms for spin gaps are not incompatible but comple- 
mentary. They result from the same phenomenon of virtual hopping between layers, which is 
of course only virtual because the direct coherent hopping at zero frequency is blocked by in- 
teraction effects. Both the spin gap phenomenon and superconductivity itself are explicable 
on the same basis. 

We may understand the superexchange phenomenon in Mott insulators by realizing that 
the orbitals of opposite-spin electrons need not be orthogonal. Thus a possible path for 
exchange of two electrons of opposite spins is for the second electron (of down spin, say) 
to hop back from site 2 to site 1, during the period while the electron of up spin, while 
nominally on site 1, has made a virtual transition to site 2. Thus if af^ and a2a- are the 
properly orthogonalized one-electron operators for the two sites, the actual one-electron 
operators are 

c+ = (a+ + ea^)/Vl + e2 
c+ = (a+ + ea+)/Vl + e2 

where 

e = t/U. 

These are not orthogonal because they need not be if they belong to opposite-spin electrons. 
The superexchange energy is then 



This picture can be renormalized, as discussed in Herring 0, by taking into account 
interactions between the two electrons along the exchange path, but the physics remains the 
same. 

J is the amplitude for interchange of spins keeping the same charge state. We may also 
ask for the pair tunneling amplitude, i.e., the amplitude for an up-spin down-spin pair to 
tunnel from one site to another. This is of course charge disfavored for atomic sites, but 
for tunneling between metallic layers there is no charge rigidity. In this case the down-spin 
electron hops in the opposite direction, and the amplitude is also 

with one hop taking place by virtue of non-orthogonality of wave functions in the two layers, 
the second by the actual frustrated one-electron matrix element. This two-particle tunneling 
process is the basis of the interlayer theory. It, like superexchange, can be described by an 
effective Hamiltonian for low energy states. 

It becomes clear that there is no sharp way to distinguish between antiferromagnetic 
superexchange and pair tunneling mechanisms for spin gaps and superconductivity if we 
look at the order parameters which can result. 

First let us set up some notation. Let be a momentum near the Fermi momentum of 
the two chains or planes, 1 and 2. Orthonormal electron operators are defined as 

= 4r = 4 (2) 

+ (1) _ + (2)+ _ + 

and these can combine into the even and odd eigen-operators of the kinetic energy, 

aj = ^ etc. (3) 
a, = etc. 

We know that the kinetic energy operator may be written 
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'Hk = ekijik^ + nk^) (4) 
+t±{al 02 + a\ a_2) 
= (efc + ti_) (a+ Oo + a_„) 
+ (efc - tx) (a;!" Oe + alg a_e) 

but in the case where interactions within the layers or chains are sufficiently strong, the 
splitting given by the last expression into even and odd eigenstates is not expressed in the 
actual state: causes primarily virtual transitions |0, as in the Mott insulator. (More 
precisely: t± causes no coherent transitions). This is an experimental fact in several of the 
cuprates, according to photoemission and infrared spectroscopy. 0J^. Instead of splitting, 
the electrons hop virtually, so that the electron operators which describe the actual eigen- 
excitations are roughly 



c+ = (a+ + ea+)/Vr + ^ (5) 



c+ = (a^ + ea^)/Vl + 

Now let us imagine that by the pair tunneling argument we have derived an order pa- 
rameter in which the electrons are paired in their separate layers: 

(O.P.)tl = {ct ct, + 4 ct^) (6) 

This order parameter is, in terms of orthonormal states, a mixture of the Strong and the 
Millis pairings; it may also be written 

(a+ ate + «o «io) 



(OP) 



IL 



+ ((a+ ate - flo (^-o)- 



Equally, we may imagine a Millis-Monien pairing 
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(c+ + 4 cU) (7) 

Again, in terms of orthogonal orbitals this is a mixture of the two pairings. In terms of 
even and odd 

(O.P.)mm = {at «-e(^^=) - {< att) (-^) (8) 

(Incidentally, neither pairing is of the "s, — s" type suggested by Scalapino and 
Yakovenko). The M-M pairing will also lead to superconductivity because of the e term 
which represents interlayer pair tunneling, if there is coupling to other layers as well as to 
other momenta, through conventional short-range interactions. 

These two pairings have the interesting feature that they reinforce each other for the 
even-even pairing, but have opposite sign for the odd-odd pairing. We suggest here that the 
effective pairing Hamiltonians coming from the two sources have approximately the same 
coefficient so that in effect we have only even-even pairing: the order parameter is 

O.P. ^ {at at,), 

since the odd-odd pairing is favored with — sign by the superexchange term and favored 
with -|- sign by the interlayer term approximately equally. 

This pairing is that appropriate to explain the selection rules observed by Keimer et al |^ 
in the scattering of neutrons by the bilayer material Y Ba2Cu^0T. The observation is that a 
pronounced peak appears, rather sharp in energy at ~ 42 mev but broader than instrumental 
resolution in momentum space, near vr, vr in the ah plane Brillouin zone. The peak appears 
only below (magnetic scattering reported above seems to have been an artifact) and 
is either simply quasiparticle pair production, or the same somewhat enhanced by excitonic 
interaction effects. More exotic explanations seem incompatible with the experimental facts, 
in particular, exotic collective excitations seem to have no reason to appear sharply below T^, 
and at an energy so close to the supposed value of 2 A. The /cjj -dependence is very compatible 
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with the idea that coherence factors forbid magnetic scattering between fc's with energy gaps 
which have the same sign, and enhance strongly those at energy gaps of opposite sign, which 
presumably (using BISCO as a model, and relying on Josephson interference measurements) 
occur near points X and F, separated by tt, tt in the zone. The observed BISCO gap would 
fit the A;||-dependence well. 

The dependence on k± in the c-direction is remarkable. This is a sinusoidal curve with a 
period given by the inverse of the interplanar spacing, showing that the scattering changes 
sign between the two planes. An equivalent statement is that scattering satisfies the selection 
rule even ^ odd as far as symmetry in the pair of planes is concerned. 

A little thought convinces one that this cannot be a coherence factor selection rule. 
That is, if the coherence factor is large between /c-points (0, vr, and 7r,0,7r), the gaps 
-^Otto = — AttOtt- But then the sign of A^qo must be opposite to either one or the other, 
leading to a second peak, which is not observed, either at tt, tt, or 0, 0, vr. We propose that 
this rule holds because at the energy gap the even state is preferentially occupied, the odd 
state empty. Then the dominant scattering process for a neutron is even odd, when the 
pair of particles created is primarily at the energy gap. 

The understanding of how this comes about requires us to go rather deeply into the BCS 
mechanism and the slight generalization that is necessary in this problem. 

In conventional BCS, the effective Hamiltonian for a single pair of states is 

^one pair = (cfc - /U) + n_fc| - 1) (9) 



+A,c+ c+,^ + A+ c 



-kl Cfct 



An irrelevant constant term has been added to make it clear that Ti. may be taken such that 
it has no effect on the subspace + n^ki = 1 of states not satisfying the Schrieffer pairing 
condition. 

In the pseudospin representation, |Ty| one may write (1) in terms of the Nambu spinors 



T as 



n = {ek- /i) + yr+ + ^T, (10) 



and it can be diagonalized by 

= Uk -^{r^ = +1) + Vk -^{r^ = -1) (11) 

with 

6 sin 9k 

Uk = COS-,Vk = 

tanb'jfc — 



efc - A* 

In the present instance, for a pair of planes we have four Fermions belonging to a given 
A;- vector, cj.j"\ c^^"'* with o; = 1, 2 a plane index. Alternatively, we may use c^^, c^^ with 

There are two calculations which may be carried out. One is in the spin-gap situation, 
where we neglect all coupling to other momenta. Then the states of momentum k, as 
a decoupled subspace, couple through the pair tunneling and superexchange interactions. 
These are, first, from (1) 



- [c 

+ (1 ^ 2) + H.C. 

and, second, as previously proposed. 



^^+(1) ^+(1) ^(2) „(2) 
2 \^k] ^-fci ^-ki ^fct 



k 



In a similar way to the BCS case, states not satisfying n — 2 are annihilated by the in- 
teraction. In this case the Hamiltonian is number- conserving overall and the state which 
diagonalizes the sum of the two interactions is simply 

^ sg f-fcl" '^-kl ^vac 

and all other states are unaffected by the pairing Hamiltonian. This is the spin gap state, 
identical to that proposed by Strong and Anderson except that the pairing is in the even 



state, while the three other N — 2 states are either repulsive or not lowered in energy. 
Otherwise the story is unchanged. In the superconducting state we presume the gap function 
is correlated by couplings to other nearby momentum states and to other planes. We can 
model this by assuming that the momentum-conserving 5-function is not exact but has a 
finite width 5 k ~ j-, (an additional length scale which enters the inter layer theory, whose 
physical consequences we do not explore here) . But L — > oo leads to physical nonsense when 
explored in too great detail. 

Note that any electron excited from a spin gap state must leave behind a hole in an even 
state c^^^ and any k in the spin gap state can only accept an electron in an odd state c^^"' ■ 
Thus the spin gap states exactly satisfy the Keimer selection rule e — > o for scattering of 
electrons by neutrons. This reflects the source of the pairing energy in the kinetic energy 
term. 



Also note that although the odd states are nominally unpaired, nonetheless adding or 
removing an electron in an odd state from a spin gapped k- value destroys the pairing criterion 
n — 2 and costs one unit of pairing energy. Thus there is a gap for all one-electron excitations, 
even though only the even state is paired. 

Now we consider the more conventional BCS-like theory. Here we must require phase- 
coherence of the pairing among all /c-states, so that the pairing Hamiltonian is no longer 
number-conserving, and may be treated by the usual mean fleld theory. But when we 
transcribe our second-order pairing Hamiltonian into "even" and "odd" language, it turns 
out to read 



Thus it is favorable for pairing to occur in even or odd states but not in both together. 
The resulting problem is more complicated than BCS; in fact, unlike BCS mean field is not 



K.E. = t±{ne - Uo). 



k k-k'<L 
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quite adequate to the solution. Basically, the system goes from 4 states empty 

(a) ={nl = nl = o) to all full (b) « = < = 2) 
via even pairing only 

(c) nl = 2,nl = o 

but continuously, with coherent occupation of (a), (b), (c). A self-consistent BCS calculation 
gives us 

\k'=k\(L 



changing sign and therefore vanishing at k = kp- JTTl What we expect will happen, however, 
is that at /c ~ /ci? there will be phase-correlated fluctuations of 6^ and b^, as in the spin gap 
state. Thus the actual energy gap does not actually vanish at any k. How to deal in any 
exact way with the phase fluctuations of the gap is as yet an unsolved problem. I speculate 
that the mathematics of gap formation may resemble more nearly that of Bose condensation 
of preformed pairs than the simple pair condensation of BCS. 

If we had only the same-layer pairing Hamiltonian which we have used in the past, i.e., 



2^ [Ck C-k' (^k' +CCJ, 



k k'-k<L 

we obtain XkxT — 1 as the effective gap equation, with Xk = this can only be satisfied 
with el {Xl- The pairing amplitude bk = Ak = sin 9k. 

cos 9k = = Uk- 

A crude approximation is to assume that this same mathematics occurs twice, once for 
even pairs and once for odd ones. Basically, goes linearly from 2 to from = — to 
0; and n% goes linearly from 2 to from e„ = to = +Afc. 

The ratio of amplitudes e — e and o — o to e ^ o can be calculated using this 
assumption naively. The ratio of occupation factors is equal to 
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I^x Kde J\ de{2 - n%) + j\ nlde j\ nl)de _ 5 
!\ n%de j\{2 - n%)de + j\ nl /,1,,„(2 - nl) 3 

This leaves out coherence factors which may considerably enhance the ratio (I estimate 
by a factor 2). This is still less skewed than the data. The state may resemble more the 
correlated "spin gap" state than this uncorrelated mean field approximation. The "spin 
gap" give a ratio of 1:0, and it is unreasonable that the superconducting state should be 
very much lower. 

In a forthcoming paper, in collaboration with S. Chakravarty, we shall show how the 
amplitude, shape and intensity of the neutron peak follows from the above ideas. Many new 
complexities are caused by the unique nature of the pairing Hamiltonian in this theory, and 
we look forward to exploring an unexpectedly rich field of physics. 

I have benefitted greatly from discussions with N.P. Ong, D.C Clark, S.P. Strong and, 
especially, S. Chakravarty. The work grew out of discussions with B. Keimer. This work 
was supported by the NSF, Grant # DMR-9104873. 
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